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ABSTRACT 

Magnetic photon splitting 7 — ► 77 , a QED process that becomes important only 
in magnetic fields approaching the quantum critical value, B cr = 4.41 x 10 13 Gauss, 
is investigated as a mechanism for attenuation of 7-rays emitted near the surface of 
strongly-magnetized pulsars. Since splitting has no threshold, it can attenuate photons 
and degrade their energies below the threshold for one-photon pair production, and in 
high enough fields it may dominate photon attenuation above pair threshold. We model 
photon splitting attenuation and subsequent splitting cascades in 7-ray pulsars, includ- 
ing the dipole field and curved spacetime geometry of the neutron star magnetosphere. 
We focus specifically on PSR1509-58, which has the highest surface magnetic field of all 
the 7-ray pulsars (Bo = 3 x 10 13 Gauss). We find that splitting will not be important for 
most 7-ray pulsars, i.e. those with Bq <J 0.2B CT , either in competition with pair produc- 
tion attenuation in pair cascades, or in photon escape cutoffs in the spectrum. Photon 
splitting will be important for 7-ray pulsars having Bq ^ 0.3B cr , where the splitting 
attenuation lengths and escape energies become comparable to or less than those for 
pair production. We compute Monte Carlo spectral models for PSR1509-58, assum- 
ing that either a full photon splitting cascade or a combination of splitting and pair 
production (depending on which splitting modes operate) attenuate a power-law input 
spectrum. We find that photon splitting, or combined splitting and pair production, 
can explain the unusually low cutoff energy (between 2 and 30 MeV) of PSR1509-58, 
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and that the model cascade spectra, which display strong polarization, are consistent 
with the observed spectral points and upper limits for polar cap emission at a range of 
magnetic colatitudes up to ~ 25°. 

Subject headings: gamma rays: pulsars: theory - pulsars: individual: (PSR1509-58) 

1. INTRODUCTION 

The discovery of at least five new 7-ray pulsars by the Compton Gamma-Ray Observatory 
(CGRO) and ROSAT has re-ignited theoretical work on the physical processes and modeling of 
high-energy radiation from pulsars. Including the previously known 7-ray pulsars, Crab (Nolan et 
al. 1993) and Vela (Kanbach et al. 1994), the recent detection of pulsed 7-rays from PSR B1706-44 
(Thompson et al. 1992), Geminga (Halpern k Holt 1992, Bertsch et al. 1992), PSR1509-58 (Wilson 
et al. 1993), PSR B1055-52 (Fierro et al. 1993), PSR B1951+32 (Ramanamurthy et al. 1995) and 
possibly PSR0656+14 (Ramanamurthy et al. 1996) bring the total to at least seven. For the first 
time, it is possible to look for the similarities and patterns in the 7-ray emission characteristics that 
may reveal clues to the origin of this emission. Examples are a preponderance of double pulses, 
a 7-ray luminosity vs. polar cap current correlation, a spectral hardness vs. characteristic age 
correlation, and spectral cutoffs above a few GeV (e.g. Thompson 1996). PSR1509-58 stands out 
among the known 7-ray pulsars as having both an unusually low spectral cutoff energy (somewhere 
between 2 and 30 MeV) and the highest inferred surface magnetic field (3 x 10 13 Gauss). It has been 
detected by the CGRO instruments operating only in the lowest energy bands, BATSE (Wilson et 
al. 1993) and OSSE (Matz et al. 1994), with the higher energy instruments, COMPTEL (Bennett 
et al. 1993) and EGRET (Nel et al. 1996), giving upper limits that require a cutoff or turnover 
between 2 and 30 MeV. There is no evidence for pulsed TeV emission (Nel et al. 1992). 

There are currently two types of models for 7-ray pulsars being investigated in detail. Polar 
cap models assume that particles are accelerated along open field lines near the neutron star by 
strong parallel electric fields (e.g. Arons 1983). The primary particles induce electromagnetic 
cascades through the creation of electron-positron pairs by either curvature radiation (Daugherty 
& Harding 1982, 1994, 1996) or inverse-Compton radiation (Sturner & Dermer 1994) 7-rays. Outer 
gap models assume that the primary particles are accelerated along open field lines in the outer 
magnetosphere, near the null charge surface, where the corotation charge changes sign, and where 
strong electric fields may develop (Cheng, Ho &: Ruderman 1986a, b; Chiang & Romani 1992; 
Romani & Yadigaroglu 1995). Since the magnetic fields in the outer gaps are too low to sustain 
one-photon pair production cascades, these models must rely on photon-photon pair production of 
7-rays, interacting with either non-thermal X-rays from the gap or thermal X-rays from the neutron 
star surface, to initiate pair cascades. 

Magnetic one-photon pair production, 7 — ► e + e _ , has so far been the only photon attenuation 
mechanism assumed to operate in polar cap cascade models. Another attenuation mechanism, 
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photon splitting 7 — > 77 , will also operate in the high-field regions near pulsar polar caps but has 
not yet been included in polar cap model calculations. The rate of photon splitting increases rapidly 
with increasing field strength (Adler 1971), so that it may even be the dominant attenuation process 
in the highest field pulsars. There are several potentially important consequences of photon splitting 
for 7-ray pulsar models. Since photon splitting has no threshold, it can attenuate photons below the 
threshold for pair production, e = 2/ sin^B, and can thus produce cutoff's in the spectrum at lower 
energies. Here 6>kB is the angle between the photon momentum and the magnetic field vectors, and e 
is (hereafter) expressed in units of m e c 2 . When the splitting rate becomes large enough, splitting can 
take place during a photon's propagation through the neutron star magnetosphere before the pair 
production threshold is crossed (i.e. before an angle ~ 2/e to the field is achieved). Consequently, 
the production of secondary electrons and positrons in pair cascades will be suppressed. Instead 
of pair cascades, one could have splitting cascades, where the high-energy photons split repeatedly 
until they escape the magnetosphere. The potential importance of photon splitting in neutron 
star applications was suggested by Adler (1971), Mitrofanov et al. (1986) and Baring (1988). Its 
attenuation and reprocessing properties have been explored in the contexts of annihilation line 
suppression in gamma-ray pulsars (Baring 1993), and spectral formation of gamma-ray bursts from 
neutron stars (Baring 1991). Photon splitting cascades have also been investigated in models of 
soft 7-ray repeaters, where they will soften the photon spectrum very efficiently with no production 
of pairs (Baring 1995, Baring & Harding 1995a, Harding & Baring 1996, Chang et al. 1996). 

In this paper we examine the importance of photon splitting in 7-ray pulsar polar cap models 
(it presumably will not operate in the low fields of outer gap models). Following a brief discussion 
of the physics of photon splitting in Section 2, we present calculations of the splitting attenuation 
lengths and escape energies in the dipole magnetic field of a neutron star. A preliminary study 
(Harding, Baring & Gonthier 1996) has shown that splitting will be the primary mode of attenuation 
of 7-rays emitted parallel to a magnetic field B ^ 0.3-B cr = 1.3 x 10 13 Gauss. We then present, 
in Section 3, photon splitting cascade models for two cases: (1) when only one mode of splitting 
(_L— >||||) allowed by the kinematic selection rules (Adler 1971, Shabad 1975) operates, suppressing 
splitting of photons of parallel polarization (so that they can only pair produce) , but still permitting 
photons of perpendicular polarization to either split once or produce pairs, and (2) when the three 
splitting modes allowed by CP (charge-parity) invariance operate, producing mode switching and a 
predominantly photon splitting cascade. In Section 4, model cascade spectra are compared to the 
observed spectrum of PSR1509-58 to determine the range of magnetic colatitude emission points (if 
any) that can produce a spectral cutoff consistent with the data. These spectra have cutoff energies 
that are decreasing functions of the magnetic colatitude. It is found that a reasonably broad range 
of polar cap sizes will accommodate the data, and that strong polarization signatures appear in 
the spectra due to the action of photon splitting. 
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2. PHOTON SPLITTING AND PAIR CREATION ATTENUATION 

The basic features of magnetic photon splitting 7 — ► 77 and the more familiar process of 
single-photon pair creation 7 — > e + e~ are outlined in the next two subsections before investigating 
their role as photon attenuation mechanisms in pulsar magnetospheres. Note that throughout 
this paper, energies will be rendered dimensionless, for simplicity, using the scaling factor m e c 2 . 
Magnetic fields will also often be scaled by the critical field B cr ; this quantity will be denoted by 
a prime: B' = B/B cr . 



The splitting of photons in two in the presence of a strong magnetic field is an exotic and com- 
paratively recent prediction of quantum electrodynamics (QED), with the first correct calculations 
of the reaction rate being performed in the early '70s (Bialynicka-Birula and Bialynicki-Birula 1970; 
Adler et al. 1970; Adler 1971). Its relative obscurity to date (compared, for example, with mag- 
netic pair creation) in the astrophysical community stems partly from the mathematical complexity 
inherent in the computation of the rate. Splitting is a third-order QED process with a triangular 
Feynman diagram. Hence, though splitting is kinematically possible, when B = it is forbidden 
by a charge conjugation symmetry of QED known as Furry's theorem (e.g. see Jauch and Rohrlich, 
1980), which states that ring diagrams that have an odd number of vertices with only external pho- 
ton lines generate interaction matrix elements that are identically zero. This symmetry is broken 
by the presence of an external field. The splitting of photons is therefore a purely quantum effect, 
and has appreciable reaction rates only when the magnetic field is at least a significant fraction of 
the quantum critical field B CT = m 2 c 3 /(eh) = 4.413 x 10 13 Gauss. Splitting into more than two 
photons is prohibited in the limit of zero dispersion because of the lack of available quantum phase 
space (Minguzzi 1961). 

The reaction rate for splitting is immensely complicated by dispersive effects (e.g. Adler 1971; 
Stoneham 1979) caused by the deviation of the refractive index from unity in the strong field. 
Consequently, manageable expressions for the rate of splitting are only possible in the limit of 
zero dispersion, and are still then complicated triple integrations (see Stoneham 1979, and also 
Baier, Mil'shtein, and Shaisultanov 1986 for electric field splitting) due to the presence of magnetic 
electron propagators in the matrix element. Hence simple expressions for the rate of splitting of 
a photon of energy w in a field B were first obtained by Bialynicka-Birula and Bialynicki-Birula 
(1970), Adler et al. (1970) and Adler (1971) in the low-energy, non-dispersive limit: uB/B CI <^ 1 . 
The total rate in this limit, averaged over photon polarizations (Papanyan and Ritus, 1972), is 
expressible in terms of an attenuation coefficient 



where a = e 2 /Tic « 1/137 is the fine structure constant, A = h/(m e c) is the Compton wavelength 



2.1. Photon Splitting Rates 
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of the electron, and #kB is the angle between the photon momentum and the magnetic field vectors. 
Here C{B') is a strong-field modification factor (derivable, for example, from Eq. 41 of Stoneham, 
1979: see Eq. below) that approximates unity when B <C B cr and scales as B~® for B 3> B cr . 

The corresponding differential spectral rate for the splitting of photons of energy uj (with 
uj <C 1 ) into photons of energies uj' and uj — uj' is 

T sp (uj,uj>) « 30 ^-^) 2 rspM . (2) 

UJ 

Equations ([l]) and (|2|) are valid (Baring, 1991) when cj-B'sin^kB 1 > which for pulsar fields and 
uj sin 6*kB ^ 2 , generally corresponds to the regime of weak vacuum dispersion. Reducing 0kB 
or B dramatically increases the photon energy required for splitting to operate in a neutron star 
environment. The produced photons emerge at an angle #kB to the field since splitting is a collinear 
process in the low-dispersion limit. 

Adler (1971) observed that in the low-energy limit, the splitting rate was strongly dependent 
on the polarization states of the initial and final photons; this feature prompted the suggestion 
by Adler et al. (1970) and Usov and Shabad (1983) that photon splitting should be a powerful 
polarizing mechanism in pulsars. The polarization-dependent rates can be taken from Eq. (23) of 
Adler (1971), which can be related to equations ([[]) or @ via 

T^sp _ 1 rpSp _ ( M\ \ rpSp _ 2M\ T sp , , 

^-MW " 2 i ll--L|l " \M 2 J ~ 3Mj + M 2 2 ' [ ' 

where the scattering amplitude coefficients 

Mi = 




are given in Adler (1971) and Eq. 41 of Stoneham (1979). In the limit of B < B cr , Mi w 26/315 
and M.2 ~ 48/315, while in the limit of B > B cv , equation (|) produces Mi « l/(6i? /3 ) and 
M2 ~ 1/(3S /4 ) . The factor of two in the numerator of the right hand side of equation (|3|) accounts 
for the duplicity of photons produced in splitting. The photon polarization labelling convention 
of Stoneham (1979) is adopted here (this standard form was not used by Adler, 1971): the label 
I refers to the state with the photon's electric field vector parallel to the plane containing the 
magnetic field and the photon's momentum vector, while _L denotes the photon's electric field 
vector being normal to this plane. Summing over the polarization modes yields the relationship for 
the strong- field modification factor in equation (|l[): 

C(B '> = ^( 2 I) 2 ( 3 - M ' + - M =) • (5) 
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Note that, in the absence of vacuum dispersion, the splitting modes _L— >_L|| , ||— >_L_L and ||— ► |||| 
are forbidden by arguments of CP (charge-parity) invariance in QED (Adler 1971); dispersive 
effects admit the possibility of non-collinear photon splitting so that there is a small but non-zero 
probability for the _L— >_L|| channel. Equations (^)— @ define the rates to be used in the analyses 
of this paper, and are valid for loB' sin #kB <C 1. The triple integral expressions that Stoneham 
(1979) derives are valid (below pair creation threshold) for a complete range (i.e. to oo ) of the 
expansion parameter cuB' sin #kB j but are not presently in a computational form suitable for use 
here. Work is in progress to address this deficiency (Baring & Harding 1996), and preliminary 
results indicate that equations ([!])— @ approximate Stoneham's (1979) formulae to better than two 
percent for tuB' sin^B < 0.2, and differs by at most a factor of around 2.5 for ujB' sinO^B ~ 1.5, 
the value relevant to the calculations of this paper; the splitting rate given by Stoneham's formulae 
initially increase above the low energy limits as luB' sinO^ increases. 

Recently there has appeared a new result on the rates of photon splitting. Mentzel, Berg 
&: Wunner (1994) presented an S-matrix calculation of the rates for the three polarization modes 
permitted by CP invariance that are considered here. While their formal development is comparable 
to an earlier S-matrix formulation of splitting in Melrose & Parle ( 1983a, b), their presentation of 
numerical results appeared to be in violent disagreement (see also their astrophysical presentation 
in Wunner, Sang & Berg 1995) with the splitting results obtained via the Schwinger proper-time 
technique by Adler (1971) and Stoneham (1979) that comprise equations ([j])-® here. These 
results have now been retracted, the disagreement being due to a sign error in their numerical 
code (Wilke & Wunner 1996). The revised results are in much better agreement with the rates 
computed by Adler (1971). However, the revised numerical splitting rates of Wilke & Wunner 
(1996) still differ by as much as a factor of 3 from Baring & Harding's (1996) computations of 
Stoneham's (1979) general formulae. Baier, Mil'shtein, &: Shaisultanov (1996) generate numerical 
results from their earlier alternative proper-time calculation (Baler, Mil'shtein, & Shaisultanov 
1986) that are in accord with Stoneham's and Adler's (1971) results and also with those of Baring 
& Harding (1996). The numerical computation of the S-matrix formalism is a formidable task. 
The proper-time analysis, though difficult, is more amenable, and has been reproduced in the limit 
of B <C B CT by numerous authors. As the S-matrix and proper-time techniques should produce 
equivalent results, and indeed have done so demonstrably in the case of pair production (see DH83 
and Tsai & Erber 1974), we choose to use the amenable proper-time results outlined above in the 
calculations of this paper. 

The above results ignore the fact that the magnetized vacuum is dispersive and birefringent, so 
that the phase velocity of light is less than c and depends on the photon polarization. Dispersion 
can therefore alter the kinematics of QED processes such as splitting (Adler 1971), and further 
dramatically complicates the formalism for the rates (Stoneham 1979). Extensive discussions of 
dispersion in a magnetized vacuum are presented by Adler (1971) and Shabad (1975); considerations 
of plasma dispersion are not relevant to the problem of gamma-ray emission from pulsars because 
they become significant only for densities in excess of around 10 27 cm~ 3 , which are only attained 



- 7- 



at the stellar surface. Adler (1971) showed that in the limit of weak vacuum dispersion (roughly 
delineated by i?'sin#kB ^ 1), where the refractive indices for the polarization states are very 
close to unity, energy and momentum could be simultaneously conserved only for the splitting 
mode -1— >|||| (of the modes permitted by CP invariance) below pair production threshold. This 
kinematic selection rule was demonstrated for subcritical fields, where the dispersion is very weak, 
a regime that generally applies to gamma-ray pulsar scenarios. Therefore, it is probable that only 
the one mode (_L— >|||| ) of splitting operates in gamma-ray pulsars. This result may be modified 
by subtle effects such as those incurred by field non-uniformity. We adopt a dual scenario in this 
paper for the sake of completeness: one in which all CP-permitted modes of splitting operate, 
and one in which Adler's kinematic selection rules are imposed. Note that in magnetar models of 
soft gamma repeaters (e.g Baring 1995, Harding and Baring 1996), where supercritical fields are 
employed, moderate vacuum dispersion arises. In such a regime, it is not clear whether Adler's 
selection rules still endure, since his analysis implicitly uses weak dispersion limits of linear vacuum 
polarization results (e.g. see Shabad 1975), and omits higher order contributions (e.g. see Melrose 
and Parle 1983a, b) to the vacuum polarization (for example, those that couple to photon absorption 
via splitting) that become prominent in supercritical fields. Furthermore, plasma dispersion effects 
may be quite pertinent to soft gamma repeater models (e.g. Bulik and Miller 1996), rendering 
them distinctly different from pulsar scenarios. 



2.2. Pair Production Rate 

One-photon pair production is a first-order QED process that is quite familiar to pulsar the- 
orists. It is forbidden in field-free regions due to the imposition of four-momentum conservation, 
but takes place in an external magnetic field, which can absorb momentum perpendicular to B. 
The rate (Toll 1952, Klepikov 1954) increases rapidly with increasing photon energy and transverse 
magnetic field strength, becoming significant for 7-rays above the threshold, co = 2/ sin#kB, and 
for fields approaching B cr . When the photon energy is near threshold, there may be only a few 
kinematically available pair states, and the rate will be resonant at each pair state threshold, pro- 
ducing a sawtooth structure (Daugherty & Harding 1983, hereafter DH83). For photon energies 
well above threshold, the number of pair states becomes large, allowing the use of a more convenient 
asymptotic expression for the polarization dependent attenuation coefficient (Klepikov 1954, Tsai 
& Erber 1974): 

1 a 

T^_ L = - J B'sm9 kB A l± (x), (6) 



HAx) « { 

where \ = {oj/2)B' sin#kB- 



' (0.31,0.15) exp(-^-) X «l 

(7) 



k (0.72, 0.48) X ~ 1/3 X»l 
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In polar cap pulsar models (e.g. Sturrock 1971, Ruderman and Sutherland 1975), high energy 
radiation is usually emitted at very small angles to the magnetic field, well below pair threshold (see 
Harding 1995, for review). The 7-ray photons will convert into pairs only after they have traveled a 
distance s comparable to the field line radius of curvature p, so that sin#kB ~ s/p. From the above 
expression, the pair production rate will be vanishingly small until the argument of the exponential 
approaches unity, i.e. when uiB' sin9^ ;> 0.2. Consequently, pair production will occur well above 
threshold when B -C 0.1i? cr and the asymptotic expression will be valid. However when B ^> 0.1-E? C r, 
pair production will occur at or near threshold, where the asymptotic expression has been shown to 
fall orders of magnitude below the exact rate (DH83). In the present calculation, we approximate 
the near-threshold reduction in the asymptotic pair production attenuation coefficient by making 
the substitution, \ — > x/E, where F = 1 + 0.42(u; sin#kB/2) -2 ' 7 m equation ([?]) (DH83). Baring 
(1988) has derived an analytic expression for the one-photon pair production rate near threshold 
which gives a result that agrees numerically with the approximation of DH83. 

Yet even the near-threshold correction to the asymptotic rate becomes poor when B ^> 0.1-B cr 
and the photons with parallel and perpendicular polarization produce pairs only (DH83) in the 
ground (0,0) and first excited (0,1) and (1,0) states respectively. Here (J, k) denotes the Landau 
level quantum numbers of the produced pairs. Therefore when the local B > 0.1-B cr , instead of 
the asymptotic form in equation (^), we use the exact, polarization-dependent, pair production 
attenuation coefficient (DH83), including only the (0,0) pair state for || polarization: 



^ PP _ a sm 



exp(-C), u;>2/sin# kB , (8) 



11 ^\Poo\ 

and only the sum of the (0,1) and (1,0) states for _L polarization: 

„ pp a sin #kB / 77 T71 111 2 \ t t\ \ 1 + (1 + 2B') 1 / 2 , m 

rf = ^ w -(E E 1 + l+p 01 )eM-0, \ [n0kB (9) 



where 



£ = (1+24) 1/2 , E l = {l+pl 1 + 2B')^ 



for 

\Pjk\ 

and 



^- sin 2 9 kB - 1 - {j + k)B' + . 
4 ' \uj sm flkB 



1/2 



e = ^7sin 2 f5 kB . (10) 



Actually, both the asymptotic and exact mean-free paths (1/T pp ) are so small in fields where 
photons pair produce at threshold that it is, in fact, not important which rate is used at very high 
field strengths (i.e. B' ;> 1). The pair production rate in this regime thus behaves like a wall at 
threshold and photons will pair produce as soon as they satisfy the kinematic restrictions on uj 
given in equations (|8|) and (^). The creation of bound pairs rather than free pairs is possible in 
fields B' ^ 0.1 (Usov h Melrose 1995), but this should not affect the present calculation since we 
do not model the full pair cascade. 
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2.3. Attenuation Lengths 

To assess the relative importance of photon splitting compared to pair production through a 
dipole magnetic field, we compute the attenuation length L , defined to be the path length over 
which the optical depth is unity: 

r(6,e) = [ L T(9 kB ,u;)ds = 1 , (11) 



o 

where ds is the pathlength differential along the photon momentum vector k and T is the attenu- 
ation coefficient for either splitting, T sp , or pair production, T pp . In this paper, attenuation lengths 
are computed as averages over polarizations of the initial photon and, for splitting, sums over the 
final polarization states. Here O^b and the photon energy k° = to are functions of the position (e.g. 
see equation [SI]), specifically measured in the local inertial frame, while 6 is the colatitude of 



emission and e is the photon energy to an observer at infinity; our treatment of curved spacetime 
is discussed immediately below. In regions where the path length is much shorter than both the 
scale length of the field strength variation or the radius of curvature of the field, L reduces to 
the inverse of the attenuation coefficient. In the calculation of the splitting attenuation lengths, 
all three CP-permitted modes are assumed to operate. The attenuation length behavior of the 
individual modes are similar. 

We assume that test photons are emitted at the neutron star surface and propagate outward, 
initially parallel or at a specified angle, #kB,o to the dipole magnetic field (see Fig. |l| for a de- 
piction of the geometry). Photon emission in polar cap models of gamma-ray pulsars can occur 
above the stellar surface (but see the discussion in Section 5), which would generate attenuation 
lengths somewhat longer than those determined here, due to the r -3 decay of the field. A surface 
origin of the photons is chosen in this paper to provide a simple and concise presentation of the 
attenuation properties. We have included the general relativistic effects of curved spacetime in a 
Schwarzschild metric, following the treatment of Gonthier <fe Harding (1994, GH94) who studied 
the effects of general relativity on photon attenuation via magnetic pair production. GH94 included 
the curved spacetime photon trajectories, the magnetic dipole field in a Schwarzschild metric and 
the gravitational redshift of the photon energy. One improvement we have made here to the treat- 
ment of GH94 is to explicitly keep track of the gravitational redshift of the photon energy as a 
function of distance from the neutron star (see Appendix for details). Our analysis is confined to 
the Schwarzschild metric because the dynamical timescales for gamma-ray pulsars are considerably 
shorter than their period (e.g. P = 0.15 sec. for PSR1509-58), so that rotation effects in the Kerr 
metric can be neglected. We have taken a neutron star mass, M = 1.4 AT© and radius, R = 10 6 cm 
in these calculations. 

Fig. H illustrates how the attenuation lengths for photon splitting and pair production vary with 
energy for different magnetic colatitudes of the emission point, for surface fields of Bq — 0.1i? cr , 
and Bq = 0.7B cr . A field of Bq = 0.7B cr is the value of the polar surface field derived from the 
magnetic dipole spin-down energy loss (Shapiro & Teukolsky 1983), using the measured P and 
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P for PSR1509-58. As noted by Usov & Melrose (1995), this is exactly twice the value of the 
surface field given by formulae in other sources (Manchester & Taylor 1977, Michel 1991), which 
assume (inaccurately) that the dipole magnetic moment /i = Bq R 3 rather than \x = Bq R 3 /2 for a 
uniformly magnetized sphere of radius R. The other 7-ray pulsars have surface field strengths in 
the range 1 — 9 x 10 12 G, or 0.02 — 0.2 B cr (the Crab and Vela pulsars have fields around 0.2 B cr ). 
Note that the attenuation lengths in Fig. ^ are for unpolarized radiation; the curves for || and _L 
polarization states look very similar. 

The curves in Fig. || have a power-law behavior at high energies, i.e. for attenuation lengths 
much less than 10 6 cm, where the dipole field is almost uniform in direction and of roughly constant 
strength. They also exhibit sharp increases at the low energy end, where photons begin to escape 
the magnetosphere without attenuation. We may estimate the behavior of the power-law portions 
of the attenuation length curves in Fig ^ as follows. Since the photons are assumed to initially 
propagate parallel to the field, the field curvature will give propagation oblique to the field only 
after significant distances are traversed, so that the obliquity of the photon to the field scales, to 
first order, as the distance travelled, sin#kB ^ s - Inserting this in equation (|l|) gives a photon 
splitting attenuation coefficient oc s 6 i.e. an optical depth tx e 5 s 7 , since T sp cx e 5 . Inversion 
then indicates that the attenuation length should vary as L oc e~ 5 / 7 : this is borne out in Fig. [2|. 
For Bq ^ 0.1-B cr , pair production occurs as soon as the threshold e t h = 2/ sin^B is crossed (cf. 
Section 2.2) during the photon propagation in the magnetosphere. Essentially, due to the enormous 
creation rate immediately above the threshold, this energy serves as an impenetrable "wall" to the 
photon. Again, since sin 9^ oc s in the early stages of propagation, the pair production attenuation 
length should scale as L oc 2/e. These proportionalities hold in both curved and flat spacetime 
since general relativistic effects distort spacetime in a smooth and differentiable manner (see the 
Appendix). However, the attenuation lengths computed in the Schwarzschild metric are about a 
factor of 1.5 lower than those computed in flat spacetime (Baring & Harding 1995b). 

The photon splitting attenuation coefficient we have used is strictly valid only below pair 
threshold. Hence, the attenuation lengths for splitting depicted in Fig. |2| can be regarded as only 
being symbolic when they exceed those for pair production, since then pair threshold is reached 
before splitting occurs. No technically amenable general expressions for the rate of splitting above 
pair threshold exist in the physics literature. But the vicinity of parameter space just below pair 
threshold is the regime of importance for 7-ray pulsar models, where the emitted photons propagate 
until they either split or they reach pair threshold, in which case they pair produce. The attenuation 
length curves near the crossover points in Fig. [2] for Bq = 0.7 B cr will require inclusion of high energy 
corrections to the attenuation coefficient (Stoneham 1979) that arise as the 7 — ► e + e~ threshold 
is approached. Currently work is in progress to compute these modifications (Baring and Harding 
1996, in preparation), and preliminary results indicate that the rate in equation (|l]) is quite accurate 
for B <^ 0.2B CT , but increases by factors of at most a few for B = 0.7B CT and uj = 2 , as mentioned 
in Section 2.1 above. 
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2.4. Escape Energies 

The energy at which the attenuation length becomes infinite defines the escape energy, below 
which the optical depth is always <C 1 , and photons escape the magnetosphere; the existence of 
such an escape energy is a consequence of the r -3 decay of the dipole field. Escape energies of 
unpolarized photons for both photon splitting and pair production are shown in Fig. |3| as a function 
of magnetic colatitude 9 of the photon emission point for different values of magnetic field strength 
(see also Harding, Baring and Gonthier, 1996). The escape energies clearly decline with 9 and are 
monotonically decreasing functions of B for the range of fields shown. The divergences as 9 — > 
are due to the divergence of the field line radius of curvature at the poles. There the maximum 
angle #kB achieved before the field falls off and inhibits attenuation is proportional to the colatitude 
6 . For photon splitting, since the rate in equation (|l|) is proportional to cu 5 sin 6 9^b , and therefore 
also the attenuation length L , it follows that the escape energy scales as e esc oc # -6//5 near the 
poles (see also Fig. ||) as is determined by the condition L ~ R . For pair production, the behaviour 
of the rate (and therefore L ) is dominated by the exponential form in equation (0) , which then 
quickly yields a dependence e esc oc 9~ 1 near the poles for Bo 0.1B cr . This behaviour extends to 
higher surface fields because production then is at threshold, which determines 

^esc OC 9 

At high fields, Bq p> 0.3B CT , there is a saturation of the photon splitting attenuation lengths and 
escape energies, due to the diminishing dependence of B in the attenuation coefficient. Likewise, 
there is a saturation of the pair production escape energy at fields above which pair production 
occurs at threshold. The pair production escape energy curves are bounded below by the pair 
threshold 2/sin^B and merge for high 9, at the pair rest mass limit, e = 2, blueshifted by the 
factor (1 - 2GM / Re 2 )' 1 / 2 ~ 1.3. Note that photon splitting can attenuate photons well below 
pair threshold. For low fields, pair production escape energies are below those for splitting, but in 
high fields, splitting escape energies are lower at all 9. The escape energies are roughly equal for 
-Bo ~ 0.3-B cr . 

The effects of curved spacetime are quite significant when compared to the attenuation lengths 
and the escape energies obtained assuming flat spacetime. A comparison of the escape energies 
for splitting and pair production, computed in flat and curved spacetime, is shown in Fig. |||. The 
largest effects are due to the increase of the surface dipole field strength by roughly a factor of 1.4, 
and the correction for the gravitational redshift of the photon, which increases the photon energy 
by roughly a factor 1.2 in the local inertial frame at the neutron star surface compared to the 
energy measured by the observer in flat space (see the Appendix). The combination of these effects 
decreases the photon splitting escape energy by a factor of about 2 compared to flat spacetime. 
The decrease in escape energy for pair production is also a factor of about 2, except at the largest 
values of 9 and B', where the pair rest mass limit is reached (cf. Fig. ||) The escape energy is then 
no longer dependent on field strength, and the ratio of the curved to flat space escape energy is 
just the redshift of the photon energy (~ 0.8) from the conversion point. This is achieved in the 
upper right hand corner of the figure; photon splitting has no such strict limit. The ratios also 
become insensitive to 9 near the poles since there the photons move almost radially, thus traveling 
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along straight trajectories, and the curved-space correction to the field is not changing rapidly with 
colatitude. The curvature of the photon trajectory in a Schwarzschild metric does not affect the 
escape energies, to first order, except in the case of emission at large colatitudes, where the photon 
wavevector makes a large angle to the radial direction. 

High energy emission from curvature radiation, inverse Compton or synchrotron by relativistic 
particles with Lorentz factor T will not beam the photons precisely along the magnetic field, but 
within some angle ~ 1/T to the field. Fig. || illustrates the effect on the escape energies of a 
non-zero angle of emission of the photons, for the case where the photons are emitted at angles 
toward the dipole axis. We have chosen the angle #kB,o = -01 rad (= 0.57°) because it would be 
the angle at which photons with e ~ 100 would be emitted through the cyclotron upscattering 
process, #kB,o — B'/e (Dermer 1990). For emission angles 6>kB,o = in Fig. |5|a, which plots e esc 
for photon splitting, e esc oc B 6 ^ 5 for Bq <^ 0.3B cr and s csc oc #~ 6 / 5 , for 8 20°, dependences that 
naturally follow from the form of equation Generally, the escape energy is insensitive to the 
emission angle for 8 ^ 106>kB,o- For small angles, the escape energy decreases and the 6>kB,o = -57° 
curves flatten below the #kB,o = curves, converging as 8^0 to an energy that is proportional 
to (Bo sin #kB,o)~ 6 ^ 5 ( see Eq. [15]). This convergence is a consequence of the field along photon 



trajectories that originate near the pole being almost uniform and tilted at about angle #kB,o to the 
photon path. In Fig. ^b, the same effect is seen for pair creation, but this time the "saturation" is at 
the redshifted threshold energy 2(1 — 2GM / Re 2 ) 1 / 2 / sin6>kB,o> an d is independent of Bq . We note 
that this behaviour at low colatitudes was observed, in the case of pair creation in flat spacetime, 
by Chang, Chen and Ho (1996). 

An obvious exception to this expected behaviour is seen in Fig. |b for the B' = 3.1 curve, where 
the escape energy is actually larger at small colatitudes ( 1° <^ 8 <J 10° ) when the emission angle 
$kB,o is increased. This counter-intuitive result can be understood with the aid of Fig. |j which 
shows the increase in sin#kBi an d ^sin#kB> determined in the local inertial frame, with path length 
s along the photon trajectory. Note that (i) the %b,o = curves increase in proportion to s when 
s/R <C 1 , as described in the Appendix, and (ii) the sin#kB curves increase logarithmically with 
s/R when s/R is not very small. In this large field, pair production occurs when the threshold 
a>sin#kB = 2 is crossed, at the same path length for both 6>kB,o = and 6*kB,o = 0.57°. The 
differences in the photon trajectories (which are almost radial) for these two cases are so small that 
s effectively represents the same height above the stellar surface for both #kB,o • Since u sin #kB ~ 2 
defines the pair creation "wall" for both photon paths, the only difference in escape energies is due 
to the factor of sin#kB at the front of the pair creation rates in equations Hence, at the 

point of pair creation, the value of sin#kB is smaller for the #kB,o = 0.57° case, and therefore the 
escape energy is larger. In flat spacetime, which is not depicted in Figs. |H| or || the crossover point 
of the sin#kB curves occurs at the same s/R value as pair threshold, so that the escape energies 
are the same at this colatitude for the two cases (this situation was also observed by Chang, Chen 
and Ho 1996). Note that as photon splitting does not have the same sudden onset as pair creation, 
it takes place over a range of path lengths, mostly around 0.1 ^ s/R 2. Over this range, 
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sin6*kB in Fig. || is generally larger for the #kB,o = 0.57° case so that the splitting escape energy is 
correspondingly shorter than for emission parallel to the field, as is evident in Fig. ||a. 

3. CASCADE SPECTRA 

Here we describe briefly our Monte Carlo simulation of photon propagation and attenuation via 
splitting and pair creation in neutron star magnetospheres, together with results for single (Section 
3.2) and multiple (Section 3.3) generations of photon splitting. 

3.1. Monte Carlo Calculation 

We model the spectrum of escaping photons from a cascade above a neutron star polar cap, 
including both photon splitting and pair production, by means of a Monte Carlo simulation. The 
free parameters specified at the start of the calculation are the magnetic colatitude 9, the angles 
8k and cj)k (see Fig. [I]), the spectrum, the height above the surface zq = r — R of the photon 
emission, and the surface magnetic field strength Bq (note that entities with subscripts '0' designate 
determination at the stellar surface). From these quantities, and assuming that = without loss 
of generality, we compute the four-vectors of the photon position and momentum that are carried 
through the computation. Injected photons are sampled from a power-law distribution, 

N(e) = N e~ a , e min <e< e max (12) 

Polarization is chosen randomly to simulate unpolarized emission; this can be altered, as desired, 
for any postulated emission mechanism. 

The path of each input photon is traced through the magnetic field, in curved spacetime, 
accumulating the survival probabilities for splitting, -P s s urv , and for pair production, P[ urv , indepen- 
dently: 

P S utv(s) = expj-r(s)} (13) 

where 

t(s) = I" T{6vs,u)ds' (14) 
Jo 

is the optical depth along the path. These survival probabilities implicitly depend on the origin ro 
of the photon and its energy e at infinity. In computing the attenuation lengths (Section 2.3), the 
photon was assumed to split when the survival probability reaches 1/e, i.e. when equation ( |IT|) is 
satisfied. In the cascade simulation, the photon may either split or pair produce. The fate of each 
cascade photon is determined as follows: if the combined survival probability, -P s s U rv-fsurv > 
where 3?i is a random number between and 1, chosen at the emission point, then the photon 
escapes; if not, then if the probability that the photon survives splitting but not pair production, 
^surv(l ~~ -fsurv)/(l ~~ ^surv-^surv) > where 3?2 is a second random number, then the photon 
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pair produces; otherwise, the photon splits. When the photon splits, the energy of one of the 
final photons is sampled from the distribution given in equation (Q) and their polarizations are 
chosen from the branching ratios given in equation (|3|). The energy of the second photon from 
the splitting is determined simply from energy conservation, since both final photons are assumed 
to be collinear in the direction of the parent photon. Each final photon is then followed in the 
same way as the injected photon, with a call to a recursive procedure that stores photon energies 
and positions through many generations of splitting. When the photon pair produces, the code 
does not follow the radiation from the pairs but simply returns to the previous cascade generation. 
For field strengths typical of gamma-ray pulsars, the pair radiation, most probably synchrotron or 
inverse Compton, will not contribute significantly at the energies near the escape energy for the 
cascades where all splitting modes operate. An exception to this may occur for supercritical surface 
fields, where synchrotron photons acquire most of the energy of their primary electrons. When all 
splitting modes operate, the number of pair production events is a small fraction of the number of 
splitting events for Bq = 0.7B cr . The cascade photons are followed through many generations of 
splitting until all of the photons either escape or pair produce. The escaping photons are binned 
in energy and polarization. 



3.2. Partial Splitting Cascade 

For pulsar applications with subcritical fields, as discussed in Section 2.1, it is probable that 
the splitting modes allowed by CP invariance are further limited by kinematic selection rules to 
only the JL— »|| || mode. This restriction may be confined to regimes of weak vacuum dispersion and 
may also depend on subtleties such as field non-uniformity. Such selection rules would effectively 
prevent splitting cascades since _L photons could split only into || photons which do not split. Here 
we compute the emergent spectra in this type of cascade, a partial cascade, where _L mode photons 
can either pair produce or split into || mode photons, while the || mode photons may only pair 
produce. There is a limit of two cascade generations: one splitting and one pair production. The 
input spectrum is a power-law (Eq. [|l2|) with the parameters: e m i n = 10~ 3 , e max = 100 and 
a = 1.6. The value of the index a is chosen to match the power-law fit of the OSSE spectrum of 
PSR1509-58 (Matz et al. 1994). The maximum energy of the input spectrum e max = 100 is chosen 
to fall above the 30 MeV maximum possible cutoff or turnover energy of the observed PSR1509-58 
spectrum. For these runs, injection of 5 - 10 million photons are required to give adequate statistics. 
The number of pairs produced relative to photons in these partial splitting cascades is obviously 
higher than in the full cascades examined in the next section. Note that in more complete gamma- 
ray pulsar models that include the pair radiation, multiple generations of splitting might still be 
possible, being interspersed with generations of conventional synchrotron/pair cascading. 

Figure ^ shows partial splitting cascade spectra in each final polarization mode for photons 
injected parallel to the local magnetic field at different magnetic colatitudes. The spectra for the two 
polarization modes are cutoff at slightly different energies, reflecting the different escape energies for 



-15- 



splitting, which cuts off the _L mode, and for pair production, which cuts off the || mode. There is 
a slight bump below the cutoff in the || mode spectrum, due to escaping photons from _L— ►|||| mode 
splitting, but only an attenuation cutoff in the perp mode spectrum. Figure || illustrates the effect 
of injecting photons at an angle (in this case #kB,o = 0.57° = 0.01 radians) to the local magnetic 
field direction, toward the magnetic dipole axis. The high-energy cutoff decreases, compared to the 
case of injection parallel to the field, only in the _L mode spectrum and not at all in the || mode 
spectrum. This behavior is due to the existence of a threshold for pair production, but not for 
splitting and can be seen from Figs. |B|a and ||b. For field strengths well above B' = 0.1, where 
photons pair produce at threshold, the pair escape energy is much less sensitive to increases in 
$kB,o than is the splitting escape energy. The partial cascade spectra therefore become more highly 
polarized at small colatitudes when #kB,o is increased. 

This effect of strong polarization, both in the energy of the spectral cutoffs and the spectral 
shape just below the cutoffs, all but disappears when photon splitting is omitted from the calcu- 
lation, thereby defining a characteristic signature of the action of 7 — > 77 . Pair production has 
much less distinctive polarization features. For example, from equations (g) and (|9j), the ratio of 
the cutoff energies at pair creation threshold between the polarization states is (1 + yl + 2B')/2 . 
For surface fields of B' = 0.7 , threshold is crossed during photon propagation in regions with much 
lower fields, typically B' ~ 0.1 , so that the spectral cutoff (or escape energy) differs only by about 
5% between polarizations; such a difference would be virtually invisible in the emission spectra. 
Clearly then, splitting is primarily responsible for polarization features shown. 

3.3. Full Splitting Cascade 

We now present model cascade spectra for the case where all three photon splitting modes 
allowed by CP invariance, _L— ►||[|, _L— >_L_L and ||— >_L||, are operating, and multiple generations of 
splitting can occur. These cascades also allow for pair production by photons of either mode. As 
noted above, for the field strength of B' = 0.7 used in the spectral models for PSR1509-58, pair 
production occurs in less than 10% of conversions. 

Figure ^ shows full splitting cascade spectra in each final polarization mode for 2 million 
photons injected parallel to the local magnetic field (in curved spacetime) at different magnetic 
colatitudes. Each cascade spectrum shows a cutoff at roughly the splitting escape energy for that 
colatitude (compare to Fig. |3|), and a bump below the cutoff from the escaping cascade photons. 
The size of the bump is a function of the number of photons attenuated above the cutoff, which 
is dependent on the ratio of the maximum input energy, e ma x, and the escape energy. For these 
models, the size of the cascade bump grows with increasing 9 because e max is held constant while 
the escape energy is decreasing. The number of splitting generations ranges from 12 when 9 = 30° 
to 3 when 9 = 2°. The size of the cascade bump at a particular 9 could of course be larger or smaller 
if £max were increased or decreased, but the positions of the cutoffs would not vary. The spectrum 
of the bump is polarized, with a well-defined zero in polarization that is a characteristic signature 
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of the splitting cascade (see Baring 1995). Note that the polarization modes have reversed their 
flux dominance in the cascade bump compared to the partial splitting cascade case. 

Although we have injected unpolarized photons in these calculations for simplicity, the relative 
flux (i.e. spectra integrated over energies) of the two polarization modes generally has a complicated 
dependence on the branching ratios for splitting defined by equation (||), due to the cascading 
process and the non-uniformity of the field. Notwithstanding, the polarization at a given energy 
does not exceed a limiting value of 3/7 (Baring 1991). The cascade spectra for injection of polarized 
photons resemble the spectra in Fig. ^, though deviations from Fig. || are not exactly proportional to 
the degree of polarization of the injection spectrum due to the inherent complexity of the interplay 
of polarization states in the cascade. 

As shown in Figure |l(J injecting photons at an angle to the local field again has a much 
larger effect at small colatitudes (i.e. for 8 <; 100#kB,o)- The high-energy cutoffs in both modes 
now decrease in energy compared to the case of injection parallel to the field, and the sizes of the 
cascade bumps are larger, both being consequences of decrease in escape energy (see Fig. ||a). This 
effect is larger at smaller colatitudes. 



4. PHOTON SPLITTING CASCADE MODELS FOR PSR1509-58 

The multiwavelength spectrum of PSR1509-58, compiled from radio to TeV energies (Thomp- 
son 1996), shows that the peak in the power output from this pulsar, as is the case for most other 



7-ray pulsars, falls in the 7-ray band. Figures 11-14 show the high energy portion of this spectrum, 



near the cutoff, which we compare with our model spectra at different emission colatitudes. No 
formal procedure for fitting the data with the model was followed, since a simple visual comparison 
demonstrating the spectral cutoff is sufficient for the scientific goals of this paper. The e 2 F(e) 
format plots equal energy per logarithmic decade and clearly demonstrates the need for a cutoff 
or sharp turnover somewhere between the highest OSSE detected point at 3 MeV and the lowest 
EGRET upper limit at 30 MeV. Although there appears to be a discontinuity between the GINGA 
data points below 100 keV and the OSSE data points, it is common for separate fits of data from 
two different detectors to produce disparate results, even in the same energy range. Furthermore, 
the e 2 F(e) format tends to magnify the differences. The difference in spectral index of the Ginga 
and OSSE fits probably indicates a true break in the power-law spectrum around 100 keV. We have 
taken the OSSE index for the input spectrum for our cascade simulation since it most accurately 
measures the observed spectrum at the energies of importance for the model. The offset between 
the Ginga and OSSE data (or their different spectral indices) does not impact the conclusions of 
this paper, since the cascade formation is determined by the photon population in the upper end of 
the OSSE range. Note that while EGRET has obtained upper limits to the pulsed emission above 
around 30 MeV, there are earlier reports of a marginal detection by COS-B (e.g.. Hartmann et 
al. 1993), with data points lying above the EGRET limits. This apparent discrepancy remains to 
be resolved, and we opt here to consider only the later and superior EGRET observations. The 
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Comptel point and limits in Figs. [Ll|-fi4| are a preliminary analysis of data from VP23 (Hermsen 
et al. 1996), showing pulsed flux at 0.75 - 1 MeV and upper limits for the pulsed interval (50%) of 
the light curve. 

The cutoffs in the model photon splitting cascade spectra in Figures PI- 14 do in fact fall in 
the energy range 3-30 MeV for colatitudes less than around 30°. At colatitudes greater than ~ 30° 
the cutoffs are lower and are in severe conflict with the OSSE data points. The standard polar 
cap half-angle in flat spacetime, sin6> = (QR/c) 1 ^ 2 , for PSR1509-58 is 2.14°. Although curved 
spacetime corrections to the magnetic dipole field tend to very slightly decrease the polar cap size 
(GH94), the polar cap may be larger than the standard size due to distortion of the field near the 
light cylinder by plasma loading (Michel 1991). The results presented here assume, for simplicity, 
a single colatitude of emission for each, i.e. a polar rim rather than an extended cap. It is easy 
to envisage that a range of polar cap emission locations will produce a convolution of the spectra 
presented here, thereby generating a spectral turnover corresponding to the maximum colatitude 
of the cap, with steeper emission extending up to a cutoff defined by the minimum colatitude. The 
EGRET upper limits cannot really discern between a sharp cutoff or a more modest turnover above 
the Comptel energy range and so cast little light on the emission as a function of colatitude when 

The partial splitting cascade spectra, shown in Figs. 11 and [l^, exhibit only modest cascade 
bumps just below the cutoff. The limits on colatitude of the model spectra are essentially deter- 
mined by the cutoff energy and are restricted by the lowest EGRET upper limit to 2° <^ 9 ^ 25° in 
the #kB,o = case, and 9 <; 25° in the #kB,o = 0.57° case, where no lower limit to the colatitude is 
imposed by the observations (see below). The model spectra for 9 = 2° and 5° are only marginally 
consistent with the upper limits. The final revision of the Comptel data for PSR1509-58 (Bennett 
et al., in preparation) may require raising the lower bounds to the colatitude of emission obtained in 
this model/data comparison. The cutoff energies of these polarization-averaged spectra are some- 
what larger than the cutoff energies of the full cascade spectra (see Figs. 13 and 14) because the 



|| mode escape energies are determined solely by pair production, whose escape energies generally 
exceed those of splitting at this field strength (see Fig. §). This is especially pronounced in the 
$kB,o = 0.57° case, due to the fact that the pair production escape energy is insensitive to the 
photon emission angle for B 0.1, as is illustrated in Fig. ||b. 



The full cascade spectra, shown in Figs. |13|and|14j, have distinctive bumps below the cutoff due 
to the redistribution of photon energies via splitting. The size of the cascade bump further limits 
the magnetic colatitudes to 9 <^ 5° to avoid conflict with the Comptel upper limits. The lowest 
EGRET upper limit restricts the colatitudes to 5° ^ 9 ^ 2° in the case of emission parallel to the 
field (Fig. |^). In the case of emission at angle 6*kB,o = 0.57° (Fig. [14]), the cutoff energy in the 



cascade spectra saturates at small 9 at an energy of 25 MeV (see Fig. ||a). Consequently there is 
no low-energy limit to 9 in this case. For larger values of #kB,o> t ne spectral cutoffs would saturate 
at larger values of 9 and at lower energies. We can estimate the dependence of this saturation 
escape energy, on #kB o an d B' using the expression for the splitting attenuation coefficient in 
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equation (|l|). Assuming that 1/T sp ~ R approximately gives the escape energy: 

£^ 0.077(5' sin # kB ,o)- 6/5 , #£0.3. (15) 

This formula quite accurately reproduces the escape energies in Figs. |13] and [L4| since they are 
only weakly dependent on R, specifically oc i? -1 / 5 . When e||* < 7.8 (i.e. 4 MeV), cascade 
spectra at all colatitudes cutoff below the lowest possible observed cutoff energy for the PSR1509- 
58 spectrum. From equation (|l5|) , this occurs, for surface emission, at #kB,o ^ 0.03 for B' = 0.7. 
Therefore, splitting cascade spectra from photons emitted at larger angles to the field will not be 
compatible with the spectrum of PSR1509-58. For emission at some distance above the surface, 
the limit on #kB,o would be higher since it depends inversely on local field strength. 



All the model spectra in Figures 11-14 assume emission at the neutron star surface. Emission 
above the surface would produce higher cutoff energies at a given colatitude, due to the decrease 
in the dipole field strength with r. The upper limits on colatitude stated above would therefore 
be less restrictive for non-surface emission. Furthermore, when the field strength at the emission 
point is B ~ 0.3-B cr (at height 30% of the neutron star radius) the splitting and pair production 
escape energies are comparable, reducing the size of the splitting cascade bumps in all cases. At 
higher altitudes above the surface, pair production dominates the photon attenuation and conven- 
tional pair cascades (e.g. Daugherty & Harding 1996) would operate. Synchrotron radiation from 
the pairs would then result in a significantly softer emergent spectrum than the input power-law 
above the cyclotron energy ( + IB' — 1 rj 280 keV at the stellar surface, lower at greater radii). 
Consequently, in order to match the observations, the input power-law would have to be harder, 
and because of the remoteness of the emission point from the stellar surface, the colatitude 9 of 
emission would have to be increased substantially. 



5. DISCUSSION 

The results of this paper demonstrate that magnetic photon splitting can have a significant 
effect on 7-ray emission from the higher field (Bq ^ 10 13 G) pulsars. It can attenuate the 7-ray 
spectrum at lower energies than magnetic pair production and will do so without the creation of 
electron-positron pairs. We have found that in low fields (Bq £ 0.3-B cr ) and #kB,o = initially, 
photon splitting attenuation lengths are never shorter than those for pair production. In high fields 
(Bq ^ 0.3.B cr ), photon splitting lengths fall below those for pair production below a certain energy 
which depends on the colatitude 6. Photon splitting escape energies fall below pair production 
escape energies for Bq ^ 0.5 B cr , so that splitting may produce an observable signature for 7-ray 
pulsars having strong magnetic fields: high energy spectral cutoffs that are quite polarization- 
dependent. While pair creation alone will also generate such cutoffs, their dependence on photon 
polarization is far diminished from when splitting is active. 

We have modeled the shape of such spectral cutoffs through simulation of photon splitting 
cascades near the neutron star surface for the case of PSR1509-58. Two types of cascades result 
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from different assumptions about the selection rules governing the photon splitting modes: the 
"full splitting cascades" occur when three modes limited only by CP selection rules operate and 
the "partial splitting cascades" occur when only one mode permitted by kinematic selection rules 
operates. In the full cascades, splitting dominates the attenuation while in the partial cascades, 
pair production ultimately limits the rate at which photon energy degrades. However, the partial 
cascades show a distinct polarization signature due to the different escape energies for splitting 
and for pair production. The resulting PSR1509-58 model spectral cutoffs due to splitting and pair 
production fall in the required range for virtually all colatitudes 25°. However, the shape of the 
spectrum of full splitting cascades, due to the large reprocessing bump, is compatible with the data 
only for a very small range of colatitudes, ^ 5°. From these results we conclude that, although 
photon splitting is capable of producing spectral cutoffs well below EGRET energies regardless of 
which selection rules govern the splitting modes, the partial splitting cascades have a much larger 
range of phase space in which to operate. 

Attenuation through magnetic pair production and photon splitting near the polar cap will 
produce 7-ray spectral cutoffs that should be roughly a function of surface magnetic field strength, 
although other parameters such as polar cap size will come into play. Thus the 7-ray pulsar 
PSR0656+14, having the second highest surface field of 9.3 x 10 12 G, should have a cutoff energy 
between that of PSR1509-58 and the other 7-ray pulsars. In fact the unusually large spectral 
index of 2.8 measured by EGRET (Ramanamurthy et al. 1996) may be a pair production/photon 
splitting cutoff. 

It is thus possible to understand why PSR1509-58, with the highest magnetic field of all the 7- 
ray pulsars, has by far the lowest spectral cutoff energy and is the only 7-ray pulsar not detected by 
EGRET. In the case of Vela (Kanbach et al. 1994), Geminga (Meyer-Hasselwander et al. 1994) and 
1055-52 (Fierro et al. 1993), the spectral cutoffs observed by EGRET at a few GeV are consistent 
with one-photon pair production cascades (Daugherty and Harding 1982, 1996). Although the 
escape energies at the neutron star surface for the spin-down fields of these pulsars (Bq ~ 2— 6 x 10 12 
G) is below 1 GeV (see Fig. ||a), curvature radiation from primary electrons at one to two stellar 
radii above the surface will have pair production escape energies of several GeV. However, when 
the surface field exceeds ~ 10 13 G, photon splitting becomes the dominant attenuation mechanism 
in the electromagnetic cascades. In addition, the primary electrons may lose energy to resonant 
Compton scattering of thermal X-rays from the neutron star surface (Sturner 1995), rather than 
to curvature radiation, limiting their acceleration to much lower energies, typically 7 ~ 100. The 
resulting upscattered 7-ray spectrum is radiated much closer to the surface and will be cut off 
by photon splitting well below the EGRET energy range. It is important to emphasize that pair 
creation acting alone suffices to inhibit GeV emission in pulsars with spin-down fields as high as 
PSR1509-58, and splitting significantly enhances the attenuation and pushes spectral cutoffs to 
lower energies. 

If resonant Compton scattering losses limit the polar cap particle acceleration energies to 
7 -C 10 6 when B ;> 10 13 G, then the primary particles will radiate 7-rays via the cyclotron 
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upscattering process or CUSP (Dermer 1990). CUSP radiation would then provide the seed photons 
for the splitting cascade. The 7-ray spectrum for this process for power-law and monoenergetic 
electrons scattering thermal blackbody X-ray photons above the neutron star surface (Daugherty 
& Harding 1989) is a power-law with maximum energy e max ~ J C B' = 2 x 10 3 B i2 /Tq (Dermer 
1990), where j c is the energy above which the electrons scatter resonantly and Tq = T/10 6 K is 
the thermal X-ray temperature. In the case of PSR1509-58 with B' = 0.7, e max ~ 10 3 /Tq. Since 
the thermal surface emission component is not observed due to the strong non-thermal spectrum 
seen at X-ray energies (Kawai 1993), Tq is not known. However, PSR1509-58 is young (~ 1000 yr) 
and probably has Tq ~ 1 — 3. We would then expect s max ~ 300 - 10 3 , compatible with our choice 
of £max — 100 for the splitting cascade models. 

A dozen or so other radio pulsars have spin-down magnetic fields above 10 13 G. These pulsars 
would, like PSR1509-58, have photon splitting dominated cascades rather than pair cascades, pro- 
ducing lower yields of electron-positron pairs. It is possible that neutron stars with extremely high 
magnetic fields, where splitting is dominant at altitudes up to several stellar radii, do not produce 
sufficient pairs for coherent radio emission, an intriguing possibility. If such neutron stars exist, 
they would constitute a new class of radio quiet, low-energy 7-ray pulsars. 
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A. CURVED SPACETIME EFFECTS 

We include here some details of our treatment of general relativistic effects on the photon 
splitting and pair production attenuation in a neutron star magnetosphere. This treatment follows 
closely that of Gonthier & Harding (1994), who examined the importance of general relativistic 
effects on one-photon pair production attenuation in a Schwarzschild metric. They found that 
several effects of curved spacetime make significant corrections to the attenuation lengths and 
escape energies for this process, namely the curvature of the photon trajectories, the redshift of 
the photon energy, and the change in the dipole magnetic field. In the present calculation, the 
first two effects are included together in the expression for the photon momentum 4-vector. The 
curved trajectory of a photon in the Schwarzschild metric is confined to a single plane, which we 
may specify as the x-z plane. For an observer at rest in the local inertial frame at a radius r in 
the Schwarzschild field, the components of the momentum 4-vector are then 

nV2 



k r 



= 

k d = (l-2n)- 1/2 e , (Al) 

where u = m/r , with m = GM/c 2 as the scaled stellar mass (i.e. the Schwarzschild radius), e is 
the photon energy as observed at infinity, 

m L 2m t * ~\ 

for a neutron star radius of R , and 8 is the initial propagation angle of the photon to the radial 
direction. The and 6 components of k are adapted from equations (13) and (14) of Gonthier & 
Harding (1994), = follows from the choice of the plane of propagation, and the f component 
can be deduced from the others using k^ = . 

At the photon emission point, we first determine the angle 5 by performing two coordinate 
transformations to put the photon momentum 3-vector in the x-z plane. The spacetime trajectory 
for that photon is then computed and stored in two tables: the first giving the value of r as a 
function of the total pathlength along the trajectory, s = cAr, from the time-of-flight in the local 
frame, 

r m / r u b du 

V* K - u2 ( l - 2u )] 1/2 (1 - 2n)V2 n 2 
which closely resembles Eq. (19) of Gonthier and Harding (1994, which instead measures time in 
the non-local observer's frame), and the second table giving the value of $ as a function of r from 
the equation, 

d ^f = ul-u\l-2u) . (A4) 



Ar = r ^ , (A3) 

^ Jm/R \ — 7/2 
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The r and # are the coordinates in the x-z plane of each point along the photon trajectory. At 
each distance s along the photon path from the emission point, we look up the value of r, and from 
r deduce the value of These values of r and i9 then define the new position and momentum 
4-vectors in the x-z plane. The inverse coordinate transformations of those described above will 
then give the position and momentum 4-vectors in the frame in which we carry out the attenuation 
length calculation. 

To describe the magnetic field in curved spacetime, we use the expression of Wasserman & 
Shapiro (1983) for the dipole field measured in the local inertial frame in a Schwarzschild metric: 
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(A5) 



Note that this expression in GH94 has a typographical error in the component. For M = 1.4M 
and R = 10 km, as used in this paper, the dipole field strength at the neutron star surface at the 
pole from equation ( |A5[) is a factor ~ 1.4 times the flat space value. 

The angle #kB> obtained by taking the dot product between the photon momentum 4- vector 
and the local dipole field, is given in curved spacetime by 



B r 

cos <9 kB = — 



1-(1 



nl/2 



+ 



B» (1 - 2u) 



1/2, 



B 



Ub 



(A6) 



Since the field components B r and B^ defined by B C urved = B r f + B^-d in equation ( |A5[) are 
differentiable functions of u = m/r , it follows from equation (^^) that #kB is also differentiable 
in r . Note also that the pathlength s defined through equation ( |A3|) gives smooth ds / dr at all 
positions above the stellar surface. Consequently, the implicit function #kB(s) has a well-behaved 
Taylor series expansion (i.e. there is no singularity) about r = R (i.e. s = 0). Hence, cases where 
we set #kB,o = result in #kB oc s along the photon trajectory for (r-JJ)/iJ<l, regardless 
of whether the spacetime is curved or flat. This proportionality is responsible for certain limiting 
behaviors in the attenuation lengths and escape energies discussed in Sections 2.3 and 2.4. 
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Fig. 1. — Schematic illustration of the neutron star dipole magnetic field geometry, used for 
determination of attenuation lengths and escape energies. The dipole field has an axis in the z- 
direction, and the photon originates at position vector ro on the neutron star surface, labelled by 
the polar angle 9 . The photon propagates in the direction of its momentum vector k that makes 
an angle #kB to the local field (#kB,o a t the emission point), and is described by polar angle 9k 
with respect to its original location ro . For all results in this paper, we arbitrarily choose a photon 
trajectory in x-z plane, corresponding to a phase <p = (see the Appendix). 
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Fig. 2. — The attenuation length L for photon splitting, defined in equation (|T]) assuming three 
CP-permitted modes, and for single photon pair production as a function of energy for photons 
initially propagating parallel to the local field (i.e. #kB,o = 0), at different colatitudes 6 on the 
neutron star surface. Two cases are depicted, namely for surface fields of (a) Bq = 0.1B cr and (b) 
Bq = 0.7B cr , the latter being the spin-down field strength for PSR1509-58. At high energies e, 
the lengths scale as e~ 5 / 7 for photon splitting and e _1 for pair production, as discussed in the 
text. The lengths are averaged over photon polarizations. 
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Fig. 3. — The escape energy (i.e. where L — > oo in Eq. pjp, in units of m e c 2 , for photon 
splitting, averaged over all modes (solid curves) and in the J_— »|| || mode only (solid dots), compared 
to the escape energy for one-photon pair production (dashed curves), both as functions of magnetic 
colatitude 9 of the emission point on the neutron star surface. These energies are obtained for 
three different surface dipole magnetic field strengths and for emission along the field (#kB,o = 0). 
The curves diverge near 8 = 0, where the field line radius of curvature becomes infinite; these 
divergences scale as #~ 6 / 5 for splitting and for pair production (see text and also Fig. [| 
below). The escape energies for each process are monotonically decreasing functions of B for 
the range of parameters shown. The escape energies are averaged over photon polarizations and 
computed using the Schwarzschild metric. 
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Fig. 4. — The ratio of escape energies in curved and flat spacetime for photon splitting and pair 
production, for emission parallel to the field as a function of magnetic colatitude 6 of the emission 
point. 
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Fig. 5. — The escape energy for (a) photon splitting and (b) pair production as a function of 
magnetic colatitude for photon emission both along B and at angle #kB,o = 0.01 radians ( = 0.57° ) 
to the field. At low magnetic colatitudes 9 <J 10#kB,o > the field curvature is so low that photon 
attenuation is insensitive to the value of 6 and is well described by the uniform field results in 
equations ([l]) and (|^i0[). The O^Bfl = (solid) curves have slopes of -6/5 (splitting) and -1 (pair 
creation) at small 9 , as discussed in the text. 
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Fig. 6. — The variation of sin #kB (left-hand scale) and u sin (right-hand scale) with path length 
s above the neutron star surface, scaled by the radius R, in curved spacetime for two different values 
°f $kB,o- 00 is t ne photon energy in the local inertial frame and the light solid horizontal line marks 
the pair threshold. Observe that O^b oc s for s/R <C 1 . The colatitude 9 = 2° and field strength 
B' = 3.1 are chosen specifically to facilitate the understanding of Fig. [El. 
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Fig. 7. — Polarized spectra for partial photon splitting cascades, assuming unpolarized power-law 
emission (of index a = 1.6) parallel to the magnetic field (#kB,o = 0), at different magnetic colati- 
tudes, 9, as labelled. Here only photons of polarization _L split, while those of either polarization 
produce pairs. The normalization of the spectrum is arbitrary. 
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Fig. 8. — Same as Fig. [7|, but for emission at angle #kB,o = 0.01 radians ( = 0.57° ) to the field. 
Spectra differ only marginally from Fig. [7[ 
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Fig. 9. — Polarized spectra for full photon splitting cascades, assuming unpolarized power-law 
emission (of index a = 1.6) parallel to the magnetic field (#kB,o = 0), at different magnetic 
colatitudes, 9. The cutoffs occur at energies comparable to the escape energies computed in Section 
2.4. Pair creation is permitted in these runs, and is generally small away from the pole. 



-35- 




Fig. 10. — Same as Fig. [| but for emission at angle #kB,o = 0.57° to the field, towards the dipole 
axis. This explores the effect of finite opening angle of emission, namely that the attenuation is 
considerably more severe than in Fig. ^ at colatitudes close to the pole. 
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Fig. 11. — Partial photon splitting cascade spectra, obtained by averaging the spectra from Fig. ^ 
over polarizations and multiplying them by e 2 , compared to the observed spectrum from PSR1509- 
58. Data points are from GINGA (Kawai, Okayasu, and Sekimoto, 1993): circles, OSSE (Matz et 
al. 1994): diamonds, COMPTEL (Hermsen et al. 1996): triangles, and EGRET (Nel et al. 1996): 
squares, and the collective display is an updated version of that in Thompson (1996). 
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Fig. 12 — 



Same as Fig. 



11 



, for the model spectra of Fig. 
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Fig. 13. — The equivalent of Fig. |Tl] for full splitting cascades, i.e. obtained by averaging the 
model spectra of Fig. |9] over polarizations and multiplying them by e 2 . 
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Fig. 14 — 



Same as Fig. 



13. 



, for the model spectra of Fig. |n]. 



